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ABSTRACT 

We compare the determinations of the angular momentum of stellar mass 
black holes via the continuum and line methods with those from diskoseismology. 
The assumption being tested is that one of the QPOs (quasi-periodic oscillations) 
in each binary X-ray source is produced by the fundamental g-mode. This should 
be the most robust and visible normal mode of oscillation of the accretion disk, 
and therefore its absence should rule out diskoseismology as the origin of QPOs. 
The comparisons are consistent with the second highest frequency QPO being 
produced by this g-mode, but are not consistent with models in which one QPO 
frequency is that of the innermost stable circular orbit. 

Subject headings: accretion, accretion disks - - black hole physics - - X-rays: 
binaries 



1. Introduction 



In principle, all adiabatic perturbations of equilibrium models of accretion disks can be 
analyzed in terms of global normal modes. The pioneering studies of Shoji Kato and his 
group and the more recent work of our group have focused on accretion disks around black 
holes, so that no complications from boundary layers are involved. In addition, the only 
QPOs with fairly stable frequencies are the high frequency ones (HFQPOs) in such sources. 
This stability implies a fixed 'cavity', which suggests gravitational trapping of the modes. 



For reviews of 'relativistic diskoseismology', see iKatol (j200lf ) and IWagonerl (120081 ). A 
short summary of some observationall y relevant results from our a nalyses of the low-lyin g 
spectrum, which consi sts of g-modes flPerez et al. Hl997h . c-modes jSilberglgit_gt~aH l200lh . 



funda mental p- modes (lOrteea- Rodriguez et al 



2002h . and other p-modes (j 



Ortega-Rodriguez et al 



20081 ) is given by lWagoner et al.l (120011 ). Local analyses (restricted radi al interval) have also 
played an important role in our understanding of these perturbations (IKato et al.ll2008l ). 
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In section 2 we summarize the foundations and results of our approach. The predictions 
of black hole spin (angular momentum) are compared with those of the 'continuum' and 
'line' methods in section 3. Section 4 contains a discussion of the viability of our model and 
conclusions. 



2. Assumptions and Predictions 



We employ a fully relativistic analysis of linear hydrodynamic perturbations of geo- 
metrically thin, optically thick accretion disks. The gravitational field of the disk and its 
perturbations are negligible compared to that of the (Kerr metric) black hole. To simplify 
the analysis, we consider barotropic disks [p = p(p)], so the buoyancy frequency vanishes. 
This should be a good approximation since the disk is turbulent, thus close to being ver- 
tically isentropic. T he equilibrium disk is taken to be described by the standard relativis- 
tic thin disk model (INoyikov fc Thornd Il973t iPage fc Thornd Il974l ) [see also Section 7 of 
Ortega- Rodriguez et al. J 2008hj. We hav e also included viscosity in the perturbation analy- 
sis (jOrtega- Rodriguez fc Wagonerl 120001 ) . The dimensionless black hole angular momentum 
parameter a = cJ/GM 2 is less than unity in absolute value, but should be positive for these 
black holes accreting material from a stellar companion. 



We have applied the general relativistic formalism that llpser fc Lindbloml (Il992l ) devel- 
oped for perturbations of purely rotating perfect fluids. They show that one can express 
the Eulerian perturbations of all physical quantities through a single function proportional 
to the pressure perturbation 5p = F(r, z) exp[i(m<j) + &t)], where a is the eigenfrequency 
and m is the axial mode number. In what follows, a < and m > 0, but the same results 
apply with both signs reversed. The corotation frequency u = a + mQ, where Q(r) is the 
angular velocity of the equilibrium disk. We have adopted a weak radial WKB approxima- 
tion (characteristic radial wavelength A r <r), except near th e corotation resonance [whic h 
lies outside the region where the relevant modes are trapped (jSilbergleit k, Wagonerl 120081 )] 
. This has turned out to be a good approximation for essentially all of the modes. 



The most robust and visible mode should be the fundamental 



m 



0, and small radial 



and vertical mode numbers) g-mode, for the following reasons. It is trapped (u> 2 < n 2 ) by 
relativistic gravity slightly below the maximum of the radial epicyclic frequency k{t), which 
is in the hottest region of the accretion disk. Its r adial extent is greater than that of the 
p-modes, as well as that of the c-mode for a > 0.1 (jWagonerl 120081 ). The inner boundary of 
the p- and c-modes is near the innermost stable circular orbit (ISCO, where k — 0), which 
implies significant leakage of these modes into the (weakly radiating) inspiral region as well 
as uncertainty of the inner boundary conditions. In addition, any nonaxisymmetric mode 
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[m ytz 0) can only be made visible by Doppler boosting or eclipses by the black hole. 



The frequency of the fundamental g-mode lies below the maximum value of k by a 
fractional amount of approximately h/r ~ c s /(rQ), where h(r) is the thickness of the disk 
and c s (r) is the speed of sound on the midplane. This splitting is proportional to the 
luminosity L since the modulated photons come from the radiation pressure dominated 
region of the disk. This is valid for 0.01 < L/L Edd < 0.5 , in which case we ob tain for the 
frequency / = \ct\/2tt the key formula [plotted in Fig. 1 of IWagoner et al.l (120011 )] 



Mf/10M Q = F(a) [1 - e(L, a)]; F — 71 A - 246 Hz (a = - 1) , e « 0.1L/L Edd , (1) 

where M is the mas s of the black hole. The weak dependence of e on a is given by 
Wagoner et al.l (120011 ); the value above is sufficiently accurate for our purposes. Initia l use 
of this relation to estimate t he spin of a few black holes was made by IWagonerl (120081 ) and 
Silbergleit k Wagoner! J2008h . 



For the l owest frequency n onaxisymmetric (m = 1) g-mode, the frequency is signifi- 
cantly higher (IPerez et al.l 119971 ): 



f(m = l)/f(m = 0) = g(a) ; g « 3.5 - 5.8 (a = - 1) 



(2) 



We also note that with increasing values of m, f(m) — > m f(ISCO), the mode location 
r m r(ISCO), and its extent Ar — > (IPerez et al.lll997l ). Thus mode leakage becomes 
important and the observed modulation decreases, so they should be much less observable. 



The orbital frequency at the ISCO is given by 

Mf(ISCO)/10M Q = H(a) ; H = 220 - 1610 Hz (a = - 1) 
This frequency plays a central role in some other models of HFQPOs. 



(3) 



The most relevant global numerical simulations of black hole accretion disks have been 
carried out by iReynolds fc Miller! (120091) with i n idea l hydrodynamics as well as magneto- 
hydrodynamics (MHD), and by lO'Neill et al.l (120091 ) within viscous hydrodynamics. The 
nonrotating black hole was represented by a modified Newtonian gravitational potential. 
The evolution was typically followed for about 10 2 orbital periods of the inner disk. From 
power spectra at many radii, the m = g-mode was seen in the hydro simulations at the 
predicted frequency and radial location. It was not seen in the MHD simulations. However, 
because of the induced MRI turbulence, it would not be expected to be seen if it was at the 
same amplitude as in the hydro simulations. Because of the limited range of <fi (with periodic 
boundary conditions) and frequency in the 3D simulations, the higher m modes could not 
have been seen. 
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In the presence of viscosity which acts hydr odynamically (via the a model) and per- 



turbatively, lOrtega-Rodrfguez &: Wagoner! ( 120001 ) found that for most (including these g-) 
modes, a viscous instability is induced, with a mode growth rate of approximately afl. For 
a ~ 0.1, it is not clear that there were enough orbits in the MHD simulations to allow 
sufficient growth. An effective turbulent viscosity (generated by the magneto-rotational in- 
stability) should be present in these (thin) accretion disks, but it is not yet known in what 
ways it acts like a hydrodynamic viscosity. 



We should note that lArras et al.l ( 120061 ) found no g-modes in their shearing-box MHD 
simulations of a limited rad ial reg i on of an accretion disk. We also note that vertically 
integrated si mulations [e. g.. IChanl (120091) ] cannot capture the g-modes. From their MHD 
simulations, iTagger fc Varnierd ( 120061 ) claim that the m = 2, 3, . . . g-modes can grow to 
dominance over those of m = 0, 1 via the Rossby wave instability. However, this requires 
a large concentration of magnetic field between the black hole and the accretion disk. In 
addition, their simulation neglected the vertical structure of the disk. There were some 
indications of the generation of p-modes within the MRI-induced turbulence, however. 



3. Comparison to Determinations of Black Hole Spin from the Continuum 

and Line Methods 



The major methods of measuring the spin of black holes have been by (a) fitting the 
continuum spectrum and (b) obtaining the shape of the gravitationally redshifted wing of an 
iron line. An example of t he use of both metho ds [which essentially determine the value of 
r isco( a )/M] is provided by lSteiner et al.l ( 1201 ll ). The continuum method requires knowledge 
of the black hole mass, the inclination angle, and the distance to the source; while the line 
method requires a knowledge of the radial dependence of the emissivity of the reflection- 
fluorescence line. Of course, both methods depend upon numerical simulations which capture 
all the relevant physical conditions in the accretion disk and 'corona'. 

In Table 1 we compare determinations of black hole spin via these two methods with that 
from the fundamental g-mode [equation (pQ)], for all sources with high-frequency QPOs and 
at least one determination of spin. We have not included the errors in the QPO frequencies, 
since they are typically smaller than the errors in M. Note that for all sources with at 
least two HFQPOs, the two highest frequ encies have a ratio very close t o 3/2, which may 
provide an important clue to their origin ( lAbramowicz &: Kluzniakil200ll ). A key source is 
GROJ1655-40, since the highest frequency (requiring a > 1) cannot be due to the g-mode. 
However, the g-mode determination of spin from the lower frequency QPO is consistent with 
those from the line method, but somewhat above that from the continuum method. Note, 
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Table 1. Determinations of Black Hole Spin 



Source 




M/Mq 


/(Hz) 


Mf/10M Q 


Spin (g-mode) 


Spin (method) 




GRS1915+ 


105 


14.0±4.4 


41 


57 ± 18 


0.00 < a < 0.27 


0.54 < a < 0.58 


(L) 








67 


94 ± 29 


0.00 < a < 0.72 


0.97 < a < 0.99 


(L) 








113 


158 ± 50 


0.51 < a < 0.98 


0.98 < a < 1.00 


(C) 








168 


235 ± 74 


0.82 < a < 1.00 






XTEJ1550- 


-564 


9.1 ±0.6 


92? 


84 ±6 


0.13 < a < 0.46 


-0.11 < a < 0.71 


(C) 








184 


167 ± 11 


0.80 < a < 0.93 


0.33 < a < 0.70 


(L) 








276 


251 ± 17 


0.98 < a < 1.00 


0.75 < a < 0.77 


(10 

v. / 


GROJ1655 


-40 


6.3 ± 0.5 


300 


189 ± 15 


0.86 < a < 0.98 


0.65 < a < 0.75 


(C) 








450 


284 ± 23 


No solution 


0.90 < a < 1.00 


(L) 














0.97 < a < 0.99 


(L) 


Cyg X-l 




14.8 ± 1.0 


135 


200 ± 14 


0.90 < a < 0.99 


0.04 < a < 0.06 


(L) 














0.92 < a < 1.00 


(C) 














0.95 < a < 0.98 


(L) 


XTEJ1650- 


-500 


5 ± 2 


250 


125 ± 50 


0.07 < a < 0.92 


0.78 < a < 0.80 


(L) 


XTEJ1859- 


f226 


> 5.4 


190 


> 103 


0.46 < a 






Hl 743-322 






165 






-0.3 < a < 0.7 


(G) 








241 











Note . - — The black hole m a sses are from (in order) 



Har 



(12011 ); iGreene et all (120011 ): lOrosz et all (I2011al. 1200411: 



aftis & Greineri (120041 ); lOrosz et al. 



Corral-Santana et al.l ( 2011 ). The 



QPO frequencies are from the compilations of iRemillard &: McClintockl ( 120061 . 1201 ll ). with 
the 3/2 ratios in boldface. The 92 Hz QPO is cont roversial. The d e termin a tions of spin 
from t he continuum (C) method are f r om ( in order) iMcClintock et al.l (120061 ); ISteiner et al. 



(|201l|); IShafee et al.l (120061): iGou et al.l feOllh : ISteiner et al. 
(L) method are f rom iBlum et al 



J2009J) (first two); 



(201 2j) , whi le those from t he lin e 



Steiner et al 



(l201lh: Miller et al.l (120091) : 



Reis et al.l (120091 ): Miller et al.l (120091 ) (next two): iFabian et al.l (120121 ): Miller et al.l (120091 ) . 
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however, that these two methods do not give compatible results. 

We therefore advance the hypothesis that the lower frequency member of the four 3/2 
pairs is produced by the g-mode, and tentatively assume that is also true for the three 
single HFQPOs. Then we see that for GRS1915+105, this assumption predicts a spin which 
is consistent with all three other determinations. For XTEJ1550-564, the g-mode spin 
determination is slightly higher than those from the other methods. For Cygnus X-l, our 
determination of a is in agreement with that from the continuum method and a recent 
re sult from the line m ethod [which also agrees with the preferred value from the line method 



of iDuro et al.l (120111 )] . The earlier result from the line method is in strong disagreement 
with the other results. For XTEJ1650-500, the g-mode result is consistent with that from 
the only other spin determination, although the value of the black hole mass is uncertain. 
Unfortunately, no mass values are available for the other source with a 3/2 pair of HFQPOs, 
H1743-322. The range of spin indicated corresponds to using a probability distribution of 
the mass, indicating 5 < M/M < 15. 

For the m = 1 nonaxisymmetric g-mode, we see from equations (CQ) and §Z§ that 
Mf/10M Q « 250 - 1430 Hz for a = - 1. We note [from equation ©] that the ISCO 
frequency has a similar range of values. Thus from the observed values in Table 1 we see 
that only very small values of a are allowed, and only for the highest frequency QPO. Thus 
we can rule out these two models from the comparison with the relatively large values of 
spin required by most of the continuum and line determinations. 



Discussion and Conclusions 



An impo rtant first s tep in including the effects of magnetic fields within diskoseismology 
was taken by iFu k Lail (120091 ) . They found (within a local WKB analysis) that a vertical 
component of magnetic field could affect the trapping of the m = g-mode when its pressure 
was at least ~ 10~ 3 times the gas pressure. The inner trapping radius of the mode became 
that of the ISCO, and the mode frequency became larger than the maximum value of n(r). 
One is then led to a scenario in which the stable QPO only exists when the magnetic field 
is small enough. However, this is consistent with the observation that the HFQP Os appear 
only sporadically jMcClintock k Remillardlbood : iRemillard k McClintockll2006h . 



A key question is whether the accretion disk can be analyzed in terms of perturbations 
(normal modes plus turbulence) to a stationary structure. The turbulence could provide the 
driving force and (anti) damping of the modes, within a weakly nonlinear analysis. One is 
naturally tempted to contemplate the outcome of the mode growth. Will enough nonlinearity 
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be induced to lead to significant mode-mode coupling and related effects? Some evidence 
that nonlinear effects should be weak is the observed lack of dependence of the HFQPO 
frequencies on their amplitude. 

What perturbations (or nonlinear structures) correspond to the QPOs not associated 
with the fundament al g-mode? One would like to e xtend analyses such as the test p article 
resona nce model of lAbramowicz fc Kluzniak! (120011) o r the wave coupling model of iKato 



( 120081 ) to this problem. This has b een done by 



couplin g between epicyclic mode s (IBlaes et al. 



iorakl (120081 ). who analyzed weak nonlinear 



20071 ) in slender tori, extending the earlier 



work of lAbramowicz et al.l (120031 ) to global normal modes. When we apply a weak coupling 
analysis to our thin disk modes, a problem arises. There is a selection rule, Yl m = over 
coupled modes. To lowest order, in which three mode couplings are considered, the only 
simple way to obtain the 3/2 ratio between the two highest frequencies is by use of the 
m g-modes. However, as shown above, su ch modes are ruled ou t observationally. We 
are exploring other possibilities. We note that iRezzolla et al.l (120031 ) obtained a 3/2 ratio 
between p-modes in small accretion tori. However, the frequencies depend upon the size of 
the torus and the speed of sound. 

We anticipate that numerical simulations will give us more insight into the physical 
conditions within the accretion disk and 'corona'. However, the fact that no simulations 
have produced the observed spectra of QPOs indicates that they are incomplete. Hopefully 
the ASTROSAT satellite (the successor to RXTE, whose archives will still be providing 
important data for analysis) will soon provide new black hole sources. In the future, the 
increased area of the proposed LOFT satellite should also help us solve some of the mysteries 
of the QPOs. More precise determinations of black hole masses and spins, as well as more 
robust determinations of QPO frequencies, will hopefully provide a more definitive test of 
diskoseismology. 

We thank Alex Silbergleit and Manuel Ortega- Rodriguez for their previous contributions 
to various aspects of this research, and the referee for supplying two useful references. 
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